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Abstract 



The Landau theory of 180° domain walls in BaTi0 3 type ferroelectric particles 
is presented . Results of exact description of domain walls in bulk enabled us to 
formulate variational approach to theory of domain walls in corresponding small 
particles. The depolarization field effects and the space-charge layers are taken 
into account in the samples of the cube form. It was found that at low temper- 
atures well known hyperbolic tangent wall profile is a good approximation for 
description of domain walls in particles. Near the transition temperature it is 
more and more appropriate to speak about two walls separating ferroelectric, 
paraelectric and ferroelectric domains correspondingly as a result of splitting 
of a single ferroelectric wall into two subwalls in small BaTiOz particles. Do- 
main wall energy density, average interwall distance and change of the dielectric 
response of thick walls in small ferroelectric particles in microcompositcs is qual- 
itatively found. Our results qualitatively describe observed dependencies better 
than those theories which exist up to date. In temperature region near transition 
from the ferroelectric to paraelectric phase in micrcomposites our quantitative 
results for response, for Curie-like transition temperature and other parameters 
may be verified. 
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1 Introduction 



The transition temperature of the ferroelectric phase - paraelectric phase phase 
transition decreases with decreasing particle diameter in BaTi0 3 [1], PbTi0 3 
[2] and in KDP [3] microcomposites. It was found experimcntaly that there is 
no phase transition bellow a critical size of a particle in the microcomposite. It 
was observed also that a higher phase transition temperature in smaller parti- 
cles sometimes occurs. In thin films of KNbC>3 [4] and TGS [5, 6] both increase 
and decrease of the transition temperature were observed with decreasing par- 
ticle size. In PbTi03 microcomposites [7] the dielectric constant decreases with 
decreasing particle size, probably due to the gradual creation of a multidomain 
state in larger particles. An unusual peak of dielectric constant versus particle 
size was observed in BaTi03 [8]. 

There are two main theoretical approaches to explanation of the size effects 
in the ferroelectric - dielectric microcomposites of the above type. The signif- 
icance of depolarization field that tries to break up the ferroelectric particle 
into domains with different polarizations. The multidomain state then exhibits 
higher permittivity, but its creation becomes difficult for small enough particles 
[9]. Another possibility is that a surface layer with different transition tem- 
perature than that in the bulk is phenomenologically considered without the 
depolarization field explicitly taken into account. Then a shift of the phase 
transition temperature in the particle is present as well as the inhomogeneous 
distribution of the polarization [10]. In [10] the static properties of such particles 
is studied by numerical calculations only. The ferroelectric phase transition in 
small spherical particles was studied in [11]. Size effects were studied analyticaly 
as concerning temperature effects and polarization profile. Our model in [11] 
did not take into account depolarization fields, assuming full compensation of 
surface charges. The dynamic susceptibility deviates from Debye-like behaviour 
then and exhibits a broadening at higher frequencies. The Curie- Weiss law is 
exhibited by the static susceptibility, there exists a divergency at the point of 
the size-driven transition. Dielectric response of microcomposite ferroelectrics 
was studied in [12]. A two-phase composite of ferroelectric and dielectric mate- 
rials was described, and studied by different effective medium theories (isolated 
particles, Maxwell-Garnet theory, the effective medium approximation and the 
Bergman representation of the effective dielectric function) . There are present 
modes due to geometric resonances which lead to new low-frequency absorp- 
tion peaks near the percolation threshold of the ferroelectric material. Those 
peaks are always located higher than the peak of the transverse polar mode 
of the bulk component. The soft mode with renormalized mode strength but 
unchanged frequency exists only in the composite with percolated ferroelectric 
clusters. Below the percolation threshold the finite clusters of the ferroelectric 
material soften only to some degree, there is no soft mode. Results of exact 
description of domain walls in bulk have been found and were preliminary pub- 
lished in [13]. They enabled us to formulate variational approach to the theory 
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of domain walls in corresponding small particles. 

In this paper we study a Landau theory of 180° domain walls in BaTiOz type 
ferroelectric particles. The depolarization field effects and the space-charge lay- 
ers are taken into account in the samples of the cube form. We have found that 
at low temperatures well known hyperbolic tangent wall profile is a good ap- 
proximation for description of domain walls in particles. Nearby the first-order 
transition temperature a free energy contribution due to paraelectric domain 
state within a wall center becomes important leading to unusually large wall 
thickness. Approaching the transition temperature it is more and more appro- 
priate to speak about two walls separating ferroelectric domains as a result of 
splitting a single ferroelectric wall into two walls with a paraelectric phase be- 
tween them. Consequently the primary ferroelectric domain wall width becomes 
very large which influences also size dependence and other characteristics of the 
ferroelectric transition to the paraelectric phase in small BaTiOz particles. Do- 
main wall energy density, average interwall distance and change of the dielectric 
response of thick walls in small ferroelectric particles is qualitatively described 
and quantitatively approximately found in this temperature region. 



2 Domain walls configurations in a BaTiO^ par- 
ticle 

The free energy density F expansion into polar mode amplitude P following [14] 
has the form 

F = -P 2 + ^P 4 + -P 6 + C | VP | 2 . (1) 

2 4 6 

where the coefficient is given by a = a (T — T 0oo ), here a is a constant, T 0oo is 
the Curie- Weiss temperature for the bulk phase transition. 

Numerical values of the coefficients of the Landau free energy expansion (1) 
of BaTiOz are given in [14]. 

Assuming the domain wall localized at x=0, the polarization P oriented in 
the z-direction and that this polarisation changes in the x direction, P(x), the 
domain- wall energy 7 is defined by 

r+00 



/+00 
AF(P(x))dx, (2) 
-00 



where AF(P(x)) = F(P(x)) — F(P) is the difference of free energy densities 
between the state with a domain- wall configuration P(x) and the equilibrium 
state with a homogeneous polarization P. 

The single domain- wall configuration P(x) in BaTiO^ was in [14] described 
approximately by the following well-known profile 

P(x) = P.tanh(^), (3) 
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with P and £ as variational parameters. The domain-wall configuration P(x) 
in a BaTiOz particle has a wall energy density which may be approximately 
calculated as in [14] by 

r+D 

7 = / AF(P(x))dx, (4) 

J-D 

with 2D a distance between two neighbouring walls. It is a variational parame- 
ter. After calculating the wall energy density 7 one can then easily incorporate 
space-charge layer effects, characterised by a parameter t, in the same way is in 
[14], too. 



3 Exact solution for wall configurations 

Let us transform the free energy part (1) describing local potential of the sixth 
order in the polarization P to a more transparent form 

V(P)=b[^--(a 2 + c^+a^}. (5) 

The coefficients a, b and c above and the free energy density expansion coeffi- 
cients a, P, 7 are connected by relations: 



•SifWI-Ti- <«) 

la V o a 



c 2 



1(3 f3 Aa 



b 



a. 



To discuss properties of the local effective potential (5) it is convenient to 
introduce a dimensionless parameter p 

P=(- a )\ (7) 

which we will assume to be a nonegative quantity. The case of negative values 
leads to a two well potential of the sixth order, which is similar to the well-known 
double-well potentials. 

The local potential V(P) has two absolute minima at ±0 and there is a local 
maximum at P=0 for p = 0. The potential V(P) has two absolute minima ±a 
and there is a local minimum at P=0, and at ±c there are two local maxima for 
< p < ^. The potential V(P) has three absolute minima at ±a and at P=0, 
in ±c there are two local maxima for p = | . The potential V(P) has only single 
absolute minimum at P=0, at ±a there two local minima and in ±c there are 
two local maxima for ^ < p < 1. For p > 1 the situation just described repeats 
with the reversed role of the a and c constants. 
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The Lagrange-Euler equation for configurations of polarisation P which ex- 
tremize the free energy (1) has the form 

2C^ = bP(P 2 -a 2 ).(P 2 -c 2 ). (8) 

We assume that a > c, the equilibrium state configurations then realizes at 
points ±a (ferroelectric phase) or at P=0 (paraelectric phase). 

Severals years ago [18] we succeeded to find exact solutions of the equation 
(8). Similar solutions in different forms were obtained also elsewhere, [19] and 
[20]. For < p < | the walls are described in the form 

Jqtanh($+R) 
P{x) = ±a V * (9) 

Jl + q - tanh 2 (f + R)) 

where R is an arbitrary constant, describing the position of the center of the 
wall, and where we define q = — g- 2 - These walls interconnect two ferroelectric 
domains with opposite polarization. 

The parameter £ describes the rate of change of the polarization at the wall 
configuration and is given by 

- 1 1 (10) 



4 180° domain wall configurations 

The exact solution (9) describes a profile of a domain wall. It is convenient to 
rewrite the form (9) into a more transparent form 

P(x)=P (x).tanh(j+R), (11) 

where the first factor in (11) we call a "modulated amplitude" Po(x) of the wall 
configuration tanh. It is given by 

Po(x) = ±a , ^ (12) 



^l + q-tanh 2 (f + R)) 



While at larger distances from the center of the wall, localized at x = for 
R=0, the "modulated amplitude" takes the equilibrium configuration values 

P (±oo) = ±o. (13) 

at the center of the wall this amplitude decreases its value to 



M X = O) = ±a-0= = Po(±oo). ] /^f p (14) 
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It is clear from (14) that the nearer we are to the phase transition the 
nearer the constant p to the value \ , and consequently the smaller " modulated 
amplitude" of the wall at the center. This decrease reflects the increasing role 
of the paraelectric phase region in the center of the wall. 

The wall thus consists of two subwalls connecting domains with nonzero 
polarization with a region with almost zero polarization. At temperatures near 
the phase transition temperature the distance d w between the two subwalls 
(walls separating the ferroelectric domains and the paraelectric phase domain 
in the center) becomes larger. One can find that at this temperature region wc 
have the distance d w given by 

= 4»<1t# <15) 

which describes logarithmic divergence of the distance d w , because the nearer 
we are to the phase transition the nearer the constant p to the value ^. 

On the other hand the lower temperature the larger the "modulated am- 
plitude" Po(x — 0), approaching the same value as the amplitude at domains 
Po(±oo). One can find that in this temperature region 

d w = i^/ip, (16) 

which is a finite quantity. 

One can easily estimate temperature region in which the condition < p < ^ 
is satisfied and thus in which very thin walls may exist the distance d w of 
which logarithmicaly divergences. Using parameters for BaTiO^ from [14] we 
found that the mentioned inequalities are equivalent with temperature region 
112.24°C < T < 121.48°C. In this region thus the ferroelectric domain wall has 
in the center the large paraelectric region and consists of two walls /subwalls/ 
connecting this region with the corresponding ferroelectric region. 

From this result it follows that at 25°C the "tanh" wall profile, used in [14] 
for approximate calculations, is an adequate approximation. At this lower tem- 
peratures this material is characterized by an effective potential with only two 
absolute minima. Near phase transition from the ferroelectric phase to the para- 
electric phase in the temperature region 112.24°C < T < 121.48°C the "tanh" 
wall profile used for approximate calculations is an inadequate approximation. 
Our exact description of the domain wall is more appropriate in this region. 

5 The dielectric response of a ferroelectric par- 
ticle 

The dielectric response of a ferroelectric particle at temperature region near the 
transition temperature, where the wall width behaves according to (16), may 
be estimated in the following way. The effect of the depolarization field and 
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of the space-charge layer may be taken into account using the same variational 
approach as in [14]. However, as the wall thickness parameter the domain 
subwall distance as a domain wall width d w instead of £ from [14] is taken. 
This approach leads to the following wall energy density 7 

The equilibrium intcrwall distance 2D between two ferroelectric walls in a 
particle of a cube form with the side width L is found to be given by D 



D = Vds^ < 18 » 

where t characterizes changes of the polarisation near the surface. We see that 
increasing temperature to the transition temperature the domain wall width 
measured by the subwall distance d w increases while the interwall distance D 
decreases. Remaining structure consists of strips of ferroelectric phase domains 
separated by wide regions (domains) of the paraelectric phase. 

If both characteristic lengths, the interwall distance characteristic quantity 
D and the subwall distance characteristic quantity d w become comparable, D w 
d w , then the ferroelectric domains become very narrow and almost vanish with 
respect to usual situation in which ferroelectric domain wall width is very small 
with respect the ferroelectric domain width. The corresponding temperature is 
slightly below the transition temperature at which p = 3. 

In the regime in which £ << D, one finds the Curie- Weiss law with renor- 
malised transition temperature 

/ n Dt , 2d w a 8c . 

a = a + 3 ' 4 ^ + [-^ + 3^- (19) 

From this expression one finds that the temperature coefficient a is multiplied 
by 1 — ^jf- giving a . 

The critical Curie- Weiss temperature becomes 

T «= T --> 4 I < 20 » 

If the wall thickness (15) increases increasing temperature there will be much 
stronger renormalization of the temperature coefficient to zero approximately 
for those wall configurations at which the the interwall distance 2D becomes 
comparable with the width parameter d w . From (20) it follows that simulta- 
neously the Curie- Weiss temperature reduces to zero quickly, the ferroelectric 
phase vanishes even before the wall thickness becomes comparable with the in- 
terwall distance characteristic length D. It is also clear that the existence of the 
ferroelectric phase becomes more sensitive to the length L of the crystal due to 
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decreased value of the temperature coefficient a and due to the domain wall 
width d w dependence of the interwall distance 2D (18). 

For thin walls (d w < -^) the transition temperature dependence coefficient 
on the particle size L is larger for the domain walls as described in our paper than 
for those described approximately by simple tanh dependence used in [14] to 
qualitatively describe transition temperature dependence on L. Experimentally 
found transition temperature dependence on L is steeper than that predicted 
by the simple tanh wall model [14]. Our exact description od domain walls 
for the bulk free energy potential in BaTiOj, describing the first order phase 
transition which is near the second order phase transition and used for the 
description of the particle domain structure as in [14] is more appropriate for real 
crystals BaTiO^. Note that for thick walls (d w > -§■) the transition temperature 
dependence coefficient on the particle size L is smaller for the domain walls 
as described in our paper than for those described approximately by simple 
tanh dependence. Thus mechanism of the increasing domain wall width due to 
splitting of the domain wall into two walls (sequence ferroelectric - paraelectric - 
ferroelectric domains) may explain much more abrupt change of the Curie- Weiss 
temperature al low L values unexplained by the spontaneous strain mechanism 
in [14]. 
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